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Abstract 

A standard result by Smale states that n dimensional strongly coop- 
erative dynamical systems can have arbitrary dynamics when restricted 
to unordered invariant hyperspaces. In this paper this result is extended 
to the case when all solutions of the strongly cooperative system are 
bounded and converge towards one of only two equilibria outside of the 
hyperplane. 

An application is given in the context of strongly cooperative systems 
of reaction diffusion equations. It is shown that such a system can have 
a continuum of spatially inhomogeneous steady states, even when all 
solutions of the underlying reaction system converge to one of only three 
equilibria. 

keywords: monotone systems, reaction diffusion systems. 

1 Introduction 



Let / : IR n — > W 1 be a C 1 vector field. A dynamical system 

^ = /i( u )> i = l---n, (1) 

is said to be strongly cooperative if the following comparison principle is satis- 
fied: whenever u(t), v(t) are two solutions such that Uj(0) < 1^(0), i — 1, . . . , n, 
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then it must hold Ui(t) < Vi(t) for every t > 0, i — 1, . . . , n. Strongly cooper- 
ative systems are canonical examples of so-called strongly monotone systems, 
which have been studied extensively by M. Hirsch, H. Smith, and others, and 
more recently by Sontag and collaborators [H [3, [11] . They have applications 
in engineering problems, as well as in the study of ecological models and gene 
regulatory networks. 

Some of the most important results about abstract strongly monotone sys- 
tems guarantee a certain behavior for the generic solution of the system. For 
instance, the well known theorem by Hirsch [61 [11] states that the generic 
bounded solution converges towards the set of equilibria E (i.e. that the set 
of initial conditions, whose solution is bounded and doesn't converge towards 
E, has measure zero). Such statements are careful to exclude a small set S of 
exceptional states, about which nothing is said. 

A simple but powerful argument originally due to Smale [10] shows the 
reason for this tendency. Smale showed that any arbitrary compactly supported 
dynamical system in R" --1 can be seen as the restriction of a certain strongly 
cooperative system in M™ to an (n — l)-dimensional hyperplane. That is, the 
dynamics of strongly cooperative systems on invariant hyperplanes can a priori 
be completely arbitrary. 

Now, Smale's argument is originally carried out for so-called competitive 
systems, and the corresponding statement for strongly cooperative systems re- 
quires considering the system in the negative time direction (see for instance 
Hirsch and Smith [7], Section 3.5). But the hyperplane in question is glob- 
ally exponentially attractive in the original example. In particular, the time- 
inverted system has unbounded solutions everywhere outside of this hyperplane. 
This can be unsatisfying in some applications, for instance in the context of 
the Hirsch generic convergence theorem which requires for its usefulness the 
boundedness of solutions. 

The main result of this paper is a theorem in the spirit of the Smale theorem, 
adapted so that the strongly cooperative embedding has bounded solutions, and 
so that every solution outside of the invariant hyperplane remains bounded 
and converges towards one of only two equilibria in R n . Define the function 
S(u) := Mi + . . . + u n and the hyperplane H := S' _1 (0). 

Theorem 1 Consider a C 2 function g 
Then there exists a C 2 function f : R n 

1. for every ueK™ and every i ^ j 

2. for every u G R: g{u) = f(u); 



: H — > H , and a compact region R C H. 
— > MJ 1 , such that 

: dfi/dujiu) > 0; 
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3. There exists P > such that every solution u{t) of (QP with S(u(0)) > 
(S(u(t)) < 0) converges towards (P, . . . , P) (towards (— P, ... — P)) as 
t — > oo; 

4- if it holds that g(u) o it < on H — R, then f has no zeros other than 
those of g and ±(P, . . . , P) . 

Item 1. is a well known sufficient condition for the strong cooperativity of 
system ([1]). An important step in the proof of this theorem is the construction 
of a 'template' strongly cooperative system with bounded solutions 

^ = M(u), (2) 

which has a continuum of equilibria along a bounded subset of H. This system 
is then suitably altered so that it equals any given function g : H — > H on the 
bounded subset of H. 



Application: Strongly Cooperative Reaction Diffusion Sys- 
tems 

As an application, consider a system of reaction diffusion equations 

Ou ■ 

—^- = d i Au i + f i (u) } * = 1 ... 7i, (3) 
at 

defined on a smooth domain Q C M m under Neumann boundary conditions, 
where di > and /j : M n — > K is a C 2 function, i = 1 . . . n. A general 
question for these systems is their relationship with the corresponding finite 
dimensional system (Q. Note that while every equilibrium of (pO) corresponds 
to an equilibrium of ([3]), there may nevertheless be equilibria of ([3]) which are 
nonhomogeneous in space, and which therefore don't naturally correspond to 
an equilibrium of ([T]). 

For instance, if n — 1 and f(u) = u(u+l)(l — u), then flTJ) has exactly three 
equilibria, but it is well known that ([3]) may have finitely many more. 

As in the finite dimensional case, the condition dfi/duj > 0, i ^ j implies a 
comparison principle for the solutions of ([3]) ( [11] , Chapter 7) . This comparison 
principle is an important tool for the analysis of reaction diffusion systems in 
a variety of cases; see for instance the book by Cantrell and Cosner [2], where 
the cooperativity condition is used in the analysis of various spatial ecological 
models. 



4 



The application of Theorem [T] addresses the existence of uncountably many 
nonhomogeneous equilibria in strongly cooperative reaction diffusion systems 
with bounded solutions. 

Theorem 2 There exists a reaction diffusion system (T3|) such that 

1. for every u G 1R" and every i ^ j : df/duj(u) > 0; 

2. every solution of (T7]j converges towards one of only three equilibria; but 

3. the set of nonhomogeneous equilibria of (T3j] has the cardinality of the 
continuum. 

The construction is carried out in the case Q = (— 7r/2,7r/2) with n = 3. 
First a non-cooperative system is constructed on a two dimensional reaction dif- 
fusion system, in which a continuum of equilibria is shown to exist (Section [3]). 
Then this system is embedded in a three dimensional cooperative system using 
Theorem [IJ 

Refer to |Tl] , Chapter 7, regarding the existence and uniqueness of solutions 
of system ([3]). We use C(fl, W 1 ) as the state space for the underlying dynamical 
system. 

2 Strongly Cooperative Embedding 

In this section we construct the Smale embedding described in the introduction. 
Our first result provides the template strongly cooperative system used in the 
proof of Theorem [IJ 

Proposition 1 There exists e > and a smooth function M : IR n — > R™ such 
that 

1. for every u EW 1 and every i ^ j: ^ i (u) > 0, 

2. for every u G H, \u\ < e: M(u) = 0, and 

3. for every solution u{t) of (U|) such that S(u(0)) > (S(u(t)) < 0), it holds 
u{t) -> (1, . . . , 1) (u(t) -> (-1, ... - 1)) as t -> oo. 

Proof. 

Consider a smooth function 9 : lR n — >• [0, 1] such that i) 8(u) = on the 
closed set {u G lR n | u G H, \u\ < 1/2}; ii) 6{u) = 1 on the closed set M"-{ti6 
R n | S(u) < 1/2, |u| < 1}; iii) < 6{u) < 1 for all other u G R n . 
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Define 7 : R — ► R to be a smooth function with the following properties: i) 
for \x\ < 1, let 7(2;) := Jx, where the constant J > will be specified below; ii) 
the zeros of 7 are exactly 0, nP, —nP for some P > 0; iii) "f(x) > — 1/ (2n), a; G 
R. It follows from ii) that 7(2;) < for x > nP, j(x) > for x < —nP. 

Define the smooth function M : R™ — ► R n : 

M i (u):=e(u)(^--u^+ 1 (S(u)), i = l,...,n. (4) 

We first show the strong monotonicity of the system ([2]). Fix i = 1 . . . n, 
and let j ^ i. For |u| < 2, |5(u)| < 1: 



dM if , 9 



0(it) 

n 



1- 7 '(W). 



But 7 / (S'(m)) = J, by assumption i) in the construction of 7. Since the 

a 

du-j 



continuous function -^-[9(u)(^^- — Ui)\ has a minimum on the compact set 



{v\ \v\ < 2, \S(v)\ < 1}, we can choose J > to be large enough that fj^-(w) > 
on this set; similarly for all other choices of 

If |w| > 2 or S(u) > 1, then 9{v) = 1 on a neighborhood of w, and it holds 
for i 7^ j: 

dMi 1 1 1 

— ^ (u ) = - +1 '{S{u ))>-- — >0, 

aitj n n In 

by assumption iii) in the construction of 7. 

In order to further understand this monotone system, write M(u) = a(u) + 
b(u), where 

Oi(x) := 6{u) ^— ^— - > := 7(S'(m)), i = 1, . . . , n. 

We calculate the dot product of these vectors: 

n 1 

o(u) o 6( u ) = 6(uWS(u)) V - Ui = 0. 

i=i 

For instance, in order to determine the zeros of M(u), it holds in particular 

that M(u) = if and only if a(u) = b(u) = 0. We have a(u) = if and only 

if either 8(u) = 0, or nui = S(u) for every i. Thus the set of zeros of a is 

(Bx/2 H H) U {(a, . . . , a) | a e R}. Also, b(u) = holds exactly on the set 
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{v | S(v) = nP,0, or —nP}. The second claim in the proposition follows from 
rescaling M by a factor of 1/P. 

Let u(t) be a solution of (j2J), and define v(t) := S(u(t))/n. It follows by the 
chain rule that 

n n 1 

m/(t) = ^ uj(t) = n-y(S(u(t)))+6(u) - s ( u )~ u i = nj(3(u(t))) = n-y(nv(t)). 

i=l i=l 72 

(5) 

Thus = 7(nv(t)). In particular, H is an unstable invariant subset for 
the system (j2J). Moreover, t>(t) remains bounded and converges towards 0, P or 
— P as £ — > oo, depending on whether ^(^(O)) = 0, S(u(0)) > 0, or £(«(())) < 
respectively 

Let now Wi(t) := Ui(t) — S(u(t))/n = Ui(t) — v(t), i = 1 . . . n (so that u>(t) 
is the projection of onto if). Then it holds 

w[{t) = u'^t) - v\t) = 6{u{t)){-S{u{t)) - Ui (t)) = -e{u(t)) Wi (t) (6) 

n 

Hence the function w(t) remains bounded for t —>■ oo, and thus u(t) = 
v(t) + w(t) is also a bounded function as t — >• oo. 

To prove the third claim, consider a solution u(t) of ^ such that ^(^(O)) > 
0. Define v(t) and as above. It holds S(u(t)) > 1/2 for all large enough 
t (since S(u(t)) = nv(t) nP > 1/2 as t — > oo). But this implies that 
6(u(t)) = 1 for large t, by the definition of 9. Hence w'(t) = —w(t) for all large 
t by ([2D, and w(t) — > as t — > oo. Since -Uj(t) = Wj(t) + f (t) for every i, it 
follows Ui(t) — > P (or «j(t) — > 1 for the rescaled system), whenever S(u(0)) > 0, 
i — 1, . . . , n. Similarly for S(u(0)) < 0. ■ 

Remark 1 Let 9{u) = 1 and 7(2:) := 5tan _1 (x(x+ 1)(1 — x)) in equation 
where 5 > is small enough that 7' (a?) > — l/(2n) on R. Then one obtains a 
simplified system which satisfies all conditions of the above result, except that 
the disc of equilibria in H collapses to the single equilibrium 0. 

Proof of Theorem CD Proof. 

We continue to use the definition of the functions 9, 7, M from the proof of 
Proposition [H For notational convenience, assume w.l.o.g. that g, R have been 
rescaled so that M(u) = on R (after the construction, the embedding / can 
be rescaled back along with g and R). 
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Let G : W n —> H be the C 2 function defined by G{u) := g(u - S{u)/n) 
Thus in particular G = g on H. 

Define the C 2 function / : K n — > IR n by the equation 



fi(u) := QMi{u) + (1 - 0(u))Gi(u), 2 = 1 



for a constant Q > which will be defined shortly. The fact that / = g on 
.R is clear since M = 0, 9{u) = 0, and g> = G on this set. 
To see that the system 

is strongly monotone, we calculate the derivative 

for fixed j ^ i. Outside of the compact support of 1 — 9, this derivative is posi- 
tive by Theorem[T]i). On the support of 1-9, both §^ and ^ r [(l-9(u))G l (u)} 
are continuous functions which attain their minimum and maximum values on 
this compact set. By Theorem [1] i) , the minimum value of ^ i is positive. Thus 

for a large enough value of Q > 0, f^-(w) > on IR n . Similarly for all other 
choices of i ^ j. 

Let u(t) be a solution of ((7|), and let v(t) := S(u(t))/n as in Theorem [H 
Then 

n n n 

nv'{t) = Qj2 M M-t)) + (l-8(<t)))J2 G Mt)) = Qj2 M i(<-t)) 

i=l i=l i=l 

= Qnj(S(u(t))) = Qnj(nv(t)), 

using the fact that S(G(u)) = on lR n and (HJ). Thus it follows that v'(t) = 
Qj(nv(t)), as in (jSJ). In particular once more, if ^(^(O)) > 0, then S(u(t))/n — > 
P as t — > oo. Since / = QM whenever S(u) > 1/2, it follows u(t) — > (P, . . . , P) 
by Theorem [TJ iii). Similarly for ^(^(O)) < 0. 

In order to address the fourth statement, we show first that M(u) o u < 0, 
for all u & H . To see this, note that for w 6 H it holds S(tt) = and thus 
Mi(u) = —9(u)ui, i = 1, . . . , n. Thus M(w) o u = —9{u)u o u = —9{u) \uf < 0. 

The fact that / has no zeros outside H other than ±(P, . . . , P) follows from 
the argument above. Suppose that u o g{u) < on H — R, and consider u G H, 
f(u) = 0. If M(u) = 0, then 9{u) = by definition of 9 and necessarily 
G(u) = g{u) = 0. If M(u) ± 0, then also 9{u) < 1, else = f(u) = M{u) ^ 0. 
In that case also u ^ R by construction of M, hence u o g{u) < 0. Therefore 



8 



f(u) ou= (QM(u) + (1 - 6{u))g{u)) o«<(l- 6{u))g{u) o u < 0, 
a contradiction. ■ 

3 Continuum of Equilibria 

In this section we construct a (non-cooperative) two dimensional reaction dif- 
fusion equation with a continuum of spatially nonhomogeneous equilibria, and 
whose corresponding reaction system ([T]) is globally attractive towards a single 
equilibrium. We begin in reverse by defining the functions which will constitute 
the nonhomogeneous equilibria, and we build the reaction function based on 
them. Consider Q = (— 7r/2,7r/2) and the function 

4>\{x) := A(cos(sina;), sin(sinx)), x E Q, 

defined on Q for every A > 0. In particular, the image of 4>\ is an arc of radius 
A spanning an angle of two radians. One can easily compute 0' A (— 7r/2) = 
0' A (7r/2) = 0. Moreover, 

(0"(x))i = A sin(x) sin(sin(x)) — A cos 2 x cos(sin x), (8) 

(4>'\(x))2 = —A sin(:r) cos(sin(x)) — A cos 2 x sin(sin x). (9) 

The vector (f)'\{x) can be thought of as the acceleration vector of <f)\(x) as 
x grows from — 7r/2 to n/2, of course, and it points towards the inside of the 
circle of radius A, except for x = ±7r/2 where it is tangential to this circle. 

Let 

A := {r(cos#,sin#) | Ai < r < A 2 , -1 < 9 < 1}, 

for fixed < Ai < A 2 in M + . Given u e A, u = r(cos 9, sin 9), we denote 
A(m) := r, x{u) := sin -1 (#). Note that these functions are well defined and 
smooth on A, and that <f)\( u )(x(u)) = u. Define the vector field a : A — >• IR 2 by 

«(«) : = -€ { uMn)). (io) 

Proposition 2 There exists a smooth function g : R 2 — > M 2 such that i) g = a 
on the set A, and ii) all solutions of the system 

nil ■ 

-^=9i(u), 2 = 1,2, (11) 
converge towards a single equilibrium. 
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Proof. 

We start with some basic facts about a(u) for u G A. It follows from 
equation ([8]) that a>i{u) > for u 2 = (since x{u) = 0). It also follows, from 
equation ([9]), that a 2 {u) > 0, a 2 {u) < 0, and a 2 (M) = whenever u 2 > 0, 
?i 2 < 0, and u 2 = 0, respectively (since x{u) > 0, a;(w) < 0, and x{u) = in 
each case). 

By ffTUl) and the definition of x(u), X(u), the function a is smooth on A. Let 
«i and a 2 be embedded into smooth functions defined on a closed neighborhood 
A' of A, in such a way that both properties in the previous paragraph still hold 
on A' for the embedding function, which we also denote by a. 

Let A" be a closed neighborhood of A'. Consider three smooth functions 
Pi, P2, P3 '■ —* [0, 1] forming a partition of unity of M 2 , in the sense that 

1. Pi(u) + p 2 {u) + p(a) = 1 for every u G M. 2 , 

2. pi(w) = 1 on A, p\{u) = on R 2 — int (A'), and pi(w) G (0, 1) otherwise. 

3. p 3 (u) = on A', p 3 (u) = 1 on M 2 - int (A"), and G (0, 1) otherwise. 

In particular, p 2 = 1 — pi on A' — A, p 2 = 1 — p 3 on A" — A'. 
Let e = (ei, 0), for a fixed e\ > such that ex > ai for every a = (ai, a 2 ) G 
A". Let # : R 2 ->• M 2 be defined by 

flf(«) := pi(u)a(w) + p 2 (u)(l, 0) + p 3 (w)(e - u). 

It is clear that g is smooth and that g e a on A by construction. We 
use the Poincare-Bendixson theorem to show that every solution of ( TlTi) must 
converge towards e. First note that g = e — u outside of a bounded subset 
of M 2 , therefore a closed circular region centered on e is invariant and attracts 
all solutions. We show below that g has only one equilibrium; therefore every 
solution must converge towards e or towards some periodic orbit. But any 
periodic orbit would contain e inside its enclosed area; see for instance [T2"] . 
Section 6.8. Moreover, g 2 (ui,0) = for every U\ G JR., since this is true for 
each of the vector fields a, u — > e — u and u — > (1,0). Therefore the U\ axis 
{(ui, 0) | Mi G M} is an invariant subset of ffTTT) ; this proves that no periodic 
solutions can exist by the principle of nonintersecting orbits [12] , and that 
every solution converges towards e. 

It remains to show that g has only one equilibrium. We show this by 
considering the various subsets of the partition in M. 2 . If u G A, it holds g(u) = 
a(u) 7^ by the comments in the beginning of the proof. If u G int(A') — A, 
then still pi(u) > 0, and for u 2 ^ 0, 
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g 2 {u) = pi(u)at 2 (u) + p 2 {u) ■ = p 1 (u)a 2 (u) ^ 0, 
by the construction of the extension of a to A'. Also, for u 2 = 0, 

gx(u) = Pi(m)«i(m) + p 2 (u) • 1 > 0. 
If u G M 2 — int (A'), such that u\ < ei, then pi(w) = 0, and 
0i (it) = p 2 («) • 1 + P3(«)(ei - Mi) > 0. 
If U\ > e, then g(u) = e — it, and g(w) 7^ unless u — e. I 

Corollary 1 Let g be as in Proposition^ and let A G [Ai,A2]. Then the 
function <p\ is an equilibrium of the system 

du ■ 

-^ = A Ui + g t {u), i = 1,2, (12) 
on C(fi,IR 2 ) ; under Neumann boundary conditions. 

Proof. By construction of a on A, it holds for every A G [Ai, A2], x G that 

«(0 A (x)) = -(j>'l{x). Thus = 0'^(x) + «(0 A (x)) = + f{(f> x (x)). It was 
calculated that A (— 7r/2) = A (— 7r/2) = 0, A > 0; this satisfies the required 

boundary conditions. I 



Proof of Theorem [2} Proof. 

We use the function g from Proposition [2] in the context of Theorem [TJ Let 
7r : IR 2 — > C M 3 be a linear, metric preserving bijection. Define a smooth 
function g : H ^ H by g := 71 o g o 7r _1 . Then use Theorem [T] to embed this 
function into a strongly cooperative system ([T]). 

Recall from the proof of Proposition [2] that g = e — u outside of a bounded 
set. If \u\ > ei, o(m) = u — e, then uog{u) = uoe — uou = U\e\ — \uf < 0. Since 
the bijection ir preserves angles, this condition is also satisfied for g outside of 
a bounded region R. Thus we can use item 4. of Theorem CD to conclude that 
/ has a unique equilibrium on H. 

It follows automatically from Theorem CD that dfi/duj > for every u and 
i 7^ j. Every solution of ([1]) outside of the H converges towards ±(P, . . . P) 
(Theorem [TJ, item 3.), and every solution of ([I]) in H converges towards ir(e) by 
Proposition [2j Thus every solution of ([1]) converges towards one of only three 
equilibria. 
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Let u{x) := (p\(x), A G [Ai, A2], be any of the nonhomogeneous equilibria of 
(fl2l) . so that u"(x) + g(u(x)) = for every x G Q. Let u(x) := 7r(u(x)). By 
evaluating tt on both sides of the previous equation we obtain 

= rcu"(x) + rcg(u(x)) = u"(x) + g(u(x)) = u"(x) + f(u(x)). 

Thus u is an equilibrium of Q. After varying A over [Ai, A2], the third statement 
follows. I 



4 Discussion 

One might imagine a strengthening of Theorem [2]) in which every solution of 
the reaction system (TTJ converges towards a single equilibrium e. A standard 
'sandwich' argument shows that this is impossible: let (pQ) be strongly cooper- 
ative and converge globally towards eel* Then all solutions of the reaction 
diffusion system (J3j) must converge towards e as well. To see this, let u(x,t) 
be any solution of ([3]), and let v < u(x,0) < w componentwise, for some 
v,w G M n . If v(t),w(t) are the solutions of ([1]) with initial conditions v,w 
respectively, then v(t),w(t) also form spatially homogeneous solutions of ([3]), 
and v(t) < u(x,t) < w(t) componentwise for all t, by cooperativity [TTJ. But 
v(t),w(t) converge towards e, therefore so does u(t,x).) 

Regarding item 3. from Theorem[2], in a paper by Kishimoto and Weinberger 
[9 J it is proved that for convex Q and given a system ([3]) such that dfi/diij > 
for i ^ j, any nonhomogeneous equilibrium must be linearly unstable. More- 
over, in a recent result by Smith, Hirsch, and the author, it is shown that under 
the same hypotheses the generic bounded solution of ([3]) converges towards a 
homogeneous equilibrium [3]. In particular, the set of nonhomogeneous equi- 
libria must also be sparse (in the sense of prevalence; see [3 [3]). This indicates 
that item 3. in Theorem [2] cannot be strenghened to a substantially larger set 
of nonhomogeneous equilibria. 

Acknowledgements the author would like to thank Hal Smith for many 
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